On the effects of surrogacy of energy dissipation 
in determining the intermittency exponent in fully developed turbulence 
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The two-point correlation function of the energy dissipation, obtained from a one-point time record 
of an atmospheric boundary layer, reveals a rigorous power-law scaling with intermittency exponent 
H ~ 0.20 over almost the entire inertial range of scales. However, for the related integral moment, 
the power-law scaling is restricted to the upper part of the inertial range only. This observation 
is explained in terms of the operational surrogacy of the construction of energy dissipation, which 
influences the behaviour of the correlation function for small separation distances. 

PACS: 47.27.Jv High-Reynolds-number turbulence; 
05. 40. -a Fluctuation phenomena, random processes, noise, and Brownian motion; 
02.50.Sk Multivariate analysis 

Time records of turbulent velocity at a single point in space, obtained using a hot-wire or a laser Doppler anemome- 
ter, are usually interpreted, via Taylor's frozen flow hypothesis, as one-dimensional spatial cuts through the flow. Ve- 
locity structure functions can be obtained readily from such observables. In addition to the velocity, other quantities 
of interest include enstrophy and energy dissipation. These quantities cannot be constructed in full from the measured 
one-point velocity time series (of one or two components of velocity) and so are replaced for further analysis by the 
so-called surrogate fields. These surrogate fields usually take the form of a single component of a many-component 
field. In this note we concentrate on the surrogacy issue of energy dissipation and discuss its impact on the extraction 
of the intermittency exponent. 

To illustrate the issue, we choose turbulence measurements in an atmospheric boundary layer, made under nominally 
steady and nearly neutral conditions, in which a hot-wire probe mounted on top of a tower recorded time-series of 
both streamwise and vertical velocity components; for details of the experimental setup, see Ref. pj. The frozen 
flow hypothesis has been applied to convert the time series into spatial cuts. Upon using the method of Ref. 0, 
the Reynolds number R\ = \J {u 2 )\/v, based on the Taylor microscale A = y/ (u 2 ) / ' {{du/dx) 2 ), was determined 
to be 9000. The angular brackets denote a temporal average throughout the paper. The estimated ratio between 
integral length L, defined through the integral of the two-point correlation of the component velocity fluctuation in 
the streamwise direction, and the dissipation scale r\ = (v 3 / (e)) 1 / 4 , is 5 x 10 4 . In units of L the record length of the 
time series is L rccor d/L = 1000. The inertial range of scales is determined by examining the scaling of the third-order 
structure function; within the inertial range so determined, the power spectra show a well-defined slope close to —5/3. 
Because of instrument and cable noise, the spectral density has some amount of noise contamination towards the 
smallest scales. In order to ensure a proper construction of the derivatives dvi/dx = {vi(x + Ax) — Vi(x))j Ax, the 
noise part has been removed from the velocity signal by using a Wiener filter. 

The true energy dissipation rate 

e{x)= V -Y,{d l v 3 +d 3 v l f, (1) 

where the indices i and j represent the coordinate axes, cannot be constructed from the recorded time series. This is 
so because only the longitudinal and transverse components v x and v y of the velocity, along the streamwise direction 
x and the normal direction y respectively, are accessible. Hence, the true energy dissipation is replaced by a surrogate 
proposal. One possibility is 

£ S urrl(z) = l5v(d x V x (x)) 2 , (2) 

which, upon assuming isotropy, equals true dissipation on the average. 

For the extraction of the intermittency exponent, the two-point correlation function based on the energy dis- 
sipation field has been proposed earlier 0. Figure H shows the normalized two-point correlator (e SUT Ti(x + 
d)£ BU T:ri{x)) / (e su „i(x)) 2 . It reveals a good scaling ~ with constant intermittency exponent, fx = 0.20, over 
the extended scaling range 15ry < d < 0.3L, covering most of the inertial range. This result is in full agreement 
with older findings on turbulent jet and atmospheric boundary layer flows Q, which have also observed a clear 
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FIG. 1: Normalized two-point correlation function of the surrogate energy dissipation field obtained via the streamwise 
velocity component in the atmosphere. The dashed line has a logarithmic slope fi = 0.2 and its extent indicates the scaling 
range. Inset: logarithmic local slope of the two-point correlation function. 



and extended scaling range behaviour for the two-point correlator with about the same value for the intermittency 
exponent. 

The two-point correlator of the surrogate energy dissipation is only one tool for the extraction of the intermittenc 
exponent. Another approach that is natural from the perspective of the multifractal picture of turbulence 0, 0, 
employs the second-order moment (ef) of the coarse-grained amplitude Si = l^ 1 J l e suir (x)dx. Figure |21 shows (ef / 
as a function of the averaging scale I. The logarithmic local slope — d\n(ef) /dlnl reveals a scale-independent scaling 
exponent, which coincides roughly with the previously found value /i = 0.20 for the intermittency exponent. However, 
a peculiarity remains: the logarithmic local slope in Fig. [21 is constant (or nearly so) only in the upper part of the 
inertial range, and not over the full range, as found for the two-point correlator. One purpose of this paper is to 
understand this puzzle. 

To proceed, we note that the second-order moment and the two-point correlator are closely related (see 0) via 

(ef) = J^dxi J^dx 2 {e sm r{xi)s sm r(x 2 )) . (3) 

The left-hand-side represents a box integral over the two-point correlator and may thus be called an integral moment. 
To see when the integral moment shows the same scaling exponent as the two-point correlator, and in what range of 
scales, we assume the simplified functional form 

c (d < 1]') 

r (x + d)e sm (x)) = { a{r,'/dy (V <d< L') (4) 
1 (d > V) 



with (e SU rr) = 1, where a(r/ = 1 guarantees continuity at the decorrelation length d = L'; the parameter c is left 
free for later purposes and rj' and L' are representative small and large length scale, respectively. Upon inserting (0J 




into © we arrive at 

" <w^70 (Sf + 2 {- Th) (T ) + " «) (*)" ■ 

valid for f]' < I < L'. The first term is the targeted scaling term. The last two terms represent corrections to rigorous 
scaling. For increasing I 3> rf they fall off faster than the scaling term. If we were to fine-tune the leading-order 
correction to zero, then c = a/(l — /i). Given that fi = 0.20, the correction to the leading order term also becomes 
very small. A closer look at Fig. ^ reveals that a constant c = <z/(l — fi) is too small to approximate the small-scale 
behaviour of the two-point correlation. Consequently, the correction terms are pronounced for small I and extend 
far into the inertial range before becoming negligible. This explains qualitatively the observed scale-dependence of 
the second-order integral moment: only in the upper part of the inertial range does the scaling term with exponent 
/j, dominate, whereas for the lower part strong deviations set in, due to the behaviour of the two-point correlation 
function for very small scales. 

Figure^reveals another apparent puzzle: as the two-point distance approaches r\ the two-point correlation increases 
stronger than is suggested by the extrapolation of scaling behaviour. This is against intuition, since the onset of 
dissipation is expected to smooth out the fine structure, instead of building it up. We now offer a tentative explanation, 
which reflects the delicate issue of the surrogacy of the energy dissipation field. The two-point correlator shown in 
Fig. is based on the surrogate energy dissipation field J2J constructed from the longitudinal velocity component. 
Other constructions, such as 

15 

£ S urr2(aO = Y^(^W y (x)) 2 , (6) 

15 

£surr3(^) = -j" [2 (d x v x (x) f + {d x v y (x) ) 2 ] , (7) 
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FIG. 3: Normalized two-point correlation function of the surrogate energy dissipation fields (5J (full line), JSJ (dotted line with 
circles), and Q (dot-dashed line). The dashed straight line has a logarithmic slope /j, = 0.2. Inset magnifies the behaviour for 
short separation distances. 

are also possible. The former is based on the transverse velocity component alone, whereas the latter combines 
longitudinal and transverse components. On average, both constructions are equal in their mean value to the true 
energy dissipation field QJ, assuming isotropy. 

Figure|31compares the two-point correlator obtained from the surrogate quantities (J2J , © and J2J ■ All three variants 
reveal a rigorous power-law scaling behaviour within the extended inertial range < d < 0.3L and the corresponding 
scaling exponents are within fi = 0.20 ± 0.01, showing little differences. Only for small two-point distances d — > rj 
the two-point correlators begin to differ. Whereas the variants based on (J2J) and JHJ) practically remain identical, the 
two-point correlations based on (JJJ) are weaker for d < IO77; see inset of Fig. In view of the simplified description 
Q this implies that the constant c is smaller for the two-point correlator based on J7J than those based on © and 
©. As a consequence the leading order corrections in the expression (JSJ for the integral moment also become smaller, 
so that the scaling term should begin to dominate even at smaller length scales I. Figure [21 confirms this view: the 
local slope of the integral moment based on the surrogate field J7J) becomes constant at smaller scales than for ; 
the upper limit on the scaling range is the same for the two cases. 

When compared to the true energy dissipation , the surrogate (0) appears to be closer to (0) than the other two 
variants @ and © . We might model the amplitude of the surrogate field in terms of the amplitude of the true field 
by the relationship 

e SUII {x)=e{x)(l + f{x)) , (8) 

where f{x) behaves as a noise with mean zero. From the defining equations © and J7J we get £ S urr3 = (^surri + 
£ SU rr2)/2 = s[l + (f\ + /2)/2], so that in comparison with /1 and the noise fluctuation / 3 = (fi + /2V2 is reduced. 
Hence, this allows us to speculate that if we add more terms from the full list of Q), the extra-strong two-point 
correlations at small separation distances d < lbrj reduce further, perhaps even vanish once the surrogate field has 
converged to the true field. A quick numerical investigation, using a shear turbulence code for R\ — 99, reveals 
that the two-point correlation functions of the surrogate quantities are identical to that of true energy dissipation, 
except for very small distances, where the surrogate field possesses extra-strong correlations. While a more detailed 
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investigation is called for, this finding indicates the importance of the subtle surrogacy issue when interpreting data. 
It appears clear that the surrogacy of the energy dissipation field restricts the rigorous scaling of the second-order 
integral moment to the upper part of the inertial range; fortunately, this leaves the rigorous scaling of the two-point 
correlation function untouched over that part of the inertial range. 

Considering the behaviour of the surrogate correlators for small separation distances in terms of the modeling 
relation the noise field amplitudes cannot be expected to be uncorrelated, i.e. (f(xi)f(x2)) (f 2 )^( x i ~ x 2)i 
but should show correlations up to some separation distance. Empirically, this seems to occur within the range 
\x\ — X2I ~ 15?7. For shorter distances, the extra-strong correlation sets in for the two-point correlators based on the 
surrogate fields. A quantity related to the noise correlations is the two-point correlation of the velocity gradient field. 
Figure ^ shows {{d x Vi{x + d)){d x Vi{x))) for the measured longitudinal (wj = v x ) and transverse (vi — v y ) velocity 
components. All two-point correlations show correlations up to d w 30?7 and become zero for larger distances. 

The main message of this simple view on the surrogacy issue of the energy dissipation field is two-fold: (i) Two- 
point correlation functions of the surrogate energy dissipation field reveal rigorous scaling that is identical to the 
two-point correlation function of the true energy dissipation for two-point distances larger than about fifteen times 
the dissipation scale, (ii) The surrogacy of the energy dissipation modifies the two-point correlations for distances 
below d » 15r) and is responsible for restricting the power-law scaling behaviour of the integral moment to the upper 
part of the inertial range. A direct consequence of statements (i) and (ii) is that for the extraction of the intermittency 
exponent it is advantageous to use the two-point correlator. In view of the extended scaling range I5?y < d < 0.3L of 
the two-point correlator of the energy dissipation, it appears that the intuitive, but phcnomenological, picture of the 
scale-invariant energy cascade contains more truth than gene rally anticipated. Without further processing such as 
Extended Self Similarity [HJ or SO (3) decomposition [Til Il2l Hflj . conventional velocity structure functions in shear 
flows do not show a rigorous power-law scaling behaviour. This assigns fields such as the energy dissipation a more 
fundamental role than that of the velocity increments. As this inference is based on data analysis alone, it does not 
offer any deep theoretical explanation but requires one. 
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